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Bosonization and density-matrix renormalization group studies of the 
Fulde-Ferrell-Larkin-Ovchinnikov phase and irrational magnetization plateaus in 

coupled chains 
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We review the properties of two coupled fermionic chains, or ladders, under a magnetic field 
parallel to the lattice plane. Results are computed by complementary analytical (bosonization) and 
numerical (density-matrix renormalization group) methods which allows a systematic comparison. 
Limiting cases such as coupled bands and coupled chains regimes are discussed. We particularly focus 
on the evolution of the superconducting correlations under increasing field and on the presence of 
irrational magnetization plateaus. We found the existence of large doping-dependent magnetization 
plateaus in the weakly-interacting and strong-coupling limits and in the non-trivial case of isotropic 
couplings. We report on the existence of extended Fulde-Ferrell-Larkin-Ovchinnikov phases within 
the isotropic t-J and Hubbard models, deduced from the evolution of different observables under 
magnetic field. Emphasis is put on the variety of superconducting order parameters present at high 
magnetic field. We have also computed the evolution of the Luttinger exponent corresponding to the 
ungapped spin mode appearing at finite magnetization. In the coupled chain regime, the possibility 
of having polarized triplet pairing under high field is predicted by bosonization. 

PACS numbers: 71.10.Pm, 74.20.Mn, 74.81.-g, 75.40.Mg 
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I. INTRODUCTION 

Low-dimensional strongly correlated systems have at- 
tracted strong interest in the last years because fluctu- 
ations and energetic competitions drive these systems 
into exotic phases. Quasi-one-dimensional or strongly 
anisotropic two-dimensional superconductors are also 
known to be good candidates for the realization of Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO) phases^ 2 *^ because 
the orbital effects of the magnetic field, which induce 
the H C 2 critical field in type-II superconductors, can be 
strongly suppressed. In a singlet superconductor without 
orbital effect, there is a competition between polarizing 
the spins of the electrons and binding them into Cooper 
pairs leading to a theoretical critical field called the Pauli 
limit or the paramagnetic limit ££££ However, this limit 
can be exceeded with an inhomogeneous order parameter 
(FFLO) which can be energetically favorable, allowing 
pockets of polarized electrons and paired electrons. Qual- 
itatively, the FFLO mechanism consists in giving singlet 
Cooper pairs a finite momentum which leads to this in- 
homogeneous superconducting order parameter. Among 
these intriguing low-dimensional systems, ladders, which 
consist of a few coupled chains, proved to display deep 
new physical behaviors and sustained considerable exper- 
imental and theoretical work. 6 

At half-filling, ladders are Mott insulators and have a 
spin gap if the number of chains is even (this gap goes to 
zero in the limit of an infinite number of coupled chains) 
and no spin gap if it is odd. The two-leg ladder has 
thus a spin-liquid ground state with exponentially de- 
caying magnetic correlations. Spin gaps can also open 
under magnetic field for rational values of the magne- 
tization per site leading to plateaus in the magnetiza- 



tion curve,' Experimental evidences of zero magnetiza- 
tion plateaus have been reported on ladder and coupled- 
dimer compound o 8 ' 9 ' 10 such as Cu2(C5rIi 2 N 2 )2Cl4, 
(C 5 Hi2N) 2 CuBr4 and (5IAP) 2 CuBr 4 -2H 2 0. Away from 
half-filling, a few systems are known to develop ir- 
rational magnetization plateaus controlled by hole 
concentration ] 11 ' 12 : 13 Furthermore, when holes are intro- 
duced into the spin-liquid two-leg ladder, these charge 
carriers generically bring the system into either a metal- 
lic or a superconducting state. The isotropic two-leg lad- 
der is known to have a wide superconducting phas o 14 ' 15 
and also a metallic phase with dominant charge-density 
wavei£ (CDW) fluctuations. Another interesting feature 
is the appearance of commensurate CDW for a com- 
mensurate hole concentration^ The theoretically pro- 
posed framework to account for superconductivity in 
doped ladders relies on magnetic fluctuations and is 
based on the resonating-valence-bond (RVB) mechanism 
for superconductivity proposed by Anderson in the con- 
text of high-T c superconductors.— Within the isotropic 
Hubbard and t-J models, singlet pairing with an un- 
conventional modified d-wave structure is found. The 
competition^ 8 , between superconductivity and CDW has 
indeed been observed in the copper oxide ladder com- 
pound Sr 14 _ :! ,Ca a ;Cu24 04i + 5 (SCCO) for which the su- 
perconducting state only appears under high pressured 
However, the mechanism responsible for superconductiv- 
ity in SCCO has not reached a full agreement yet. In this 
context, the upper critical magnetic fields determined 
from transport measurement a 20 ' 21 suggest that the Pauli 
limit is exceeded in SCCO which reassesses the issue of 
the nature of the pairing. Note that superconductivity 
has been discovered in the zig-zag ladder subsystem of 
the copper oxide compound Pr2Ba4Cu70i5_a at ambi- 
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ent pressured 

Recently, a superconducting two-leg t- J ladder under a 
strong magnetic field in the plane of the ladder was stud- 
ied numerically 23 using the density-matrix renormaliza- 
tion group 24 ' 25 (DMRG) method. The magnetic curve 
displays a doping-dependent magnetization plateau, as 
predicted 2 ^ by Cabra el al. for a Hubbard ladder. In 
addition to this non trivial magnetic behavior, an ex- 
ceeding of the Pauli limit was found. Within the t-J 
model, this exceeding of the Pauli limit was explained by 
a one-dimensional analogue of the FFLO phase, hence 
reconciling the expectation of singlet pairing and the ex- 
ceeding of the Pauli limit. Lastly, the behavior of the 
superconducting correlations studied in Ref. [23[ showed 
a surprising behavior in and outside the plateau with a 
notable emergence of S z — triplet superconducting cor- 
relations. 

In this paper, we propose to greatly complete Ref. [23| 
by an extensive comparison of bosonization and numer- 
ical calculations and by extending the results to discuss 
the case of weakly-interacting or strongly anisotropic 
ladders and also to the case of weakly-coupled chains. 
In particular, we show that the picture developed in 
Ref. [23[ is consistent with bosonization and extends to 
both Hubbard and t-J models, supporting its general- 
ity. While most studies on the FFLO phase resort to 
mean-field theories for low-dimensional and unconven- 
tional superconductors^ we here use approaches which 
are more relevant for quasi one-dimensional strongly- 
correlated systems. In the same spirit, an early study on 
the superconducting phase of the t-J chain which gave ev- 
idence of developing FFLO-like correlations 2 ^ was based 
on exact diagonalization computations, and bosonization 
has also been used in this context on an attractive Hub- 
bard system. 

The article is organized as follow: first, after intro- 
ducing microscopic models (Sec. |H]), we briefly examine 
the non interacting and strong-coupling limits in Sec. IIIII 
where large doping-dependent magnetization plateaus 
can occur. The coupled bands regime is discussed in 
detail in Sec. IIVI and corresponds to the case of doped 
isotropic ladders which is the most studied at zero mag- 
netic field. Lastly, the coupled chains regime is studied 
under magnetic field (Sec. N} . 



II. MICROSCOPIC MODELS AND 
CONVENTION 

We describe the ladder system with a standard one- 
orbital Hubbard model which can have different hopping 
terms along the chains (||) and between the chains (_L). 
We consider a situation where no magnetic orbital effect 
is present and thus only keep a Zeeman coupling to the 
spin degree of freedom. This would experimentally cor- 
respond to the situation where the magnetic field H is in 
the plane of the ladder (and even along the direction of 
the ladder to minimize all orbital effects) . Then, one can 




FIG. 1: Splitting of the bands dispersion in a non-interacting 
doped ladder (/i is the chemical potential) due to Zeeman 
effect at low magnetic field H. 

write the Hubbard Hamiltonian 

i,p—l,2 : (T 
i,p=l,2 

- E H ' S ^ (!) 

i,p=l,2 

where c ip (7 , Si lP and rii lPtCr are respectively electron cre- 
ation, spin and density operators at site i on chain p and 
a is the spin index. The q\xb prefactor has been ab- 
sorbed in the definition of H for convenience. Free bands 
dispersion is given by 

£k y ,a(k) = -2t\\ cos(k) - t± cos(fcy) - H^ (2) 

(with k = k x for simplicity and k y — 0, it) and are 
sketched in FIG. [T] If we denote by k F k the Fermi wave- 
vectors, we have the Luttinger sum rule 

(ky ,a) occ. 

where n — N e / (2L) = 1 — 6 is the electron density and <5 
the hole density and L the length of the ladder. Experi- 
mentally, copper-oxide systems have a fixed hole doping 5 
rather than a fixed chemical potential fi so S will be kept 
fixed in this paper. Similarly, the magnetization per site 
m = (iV T - N l )/(2L) satisfies 

(k y ,cr) occ. 

Relations ([3]) and (|4]) are also valid in the presence of 
interactions. For quasi one-dimensional interacting sys- 
tems, a useful theorem is the generalization of the Lieb- 
Schultz-Mattis theorem to doped and magnetized states 
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done by Yamanaka-Oshikawa- Affleck 29 (YOA). This the- 
orem relates gap openings to commensurability condi- 
tions. When dealing with spinful fermions, the demon- 
stration is based on the definition of 2 twist operators 
for spins j, J, which gives two commensurability condi- 
tions for parameters in each sector. In the case of doped 
two-leg ladders, YOA parameters take the simple form2S 
1 — S + am. If a parameter is non-integer, the spectrum 
has gapless excitations in this sector. If a parameter is 
non-integer but rational, a gap can open in the sector (de- 
pending on the relevance of interactions) but necessarily 
the ground-state is degenerate with translational sym- 
metry breaking. When the parameter is an integer, the 
sector is gapless or it is gapped but with a non-degenerate 
ground state. It is interesting to note that for m, S =/= 0, 
the sector f can get gapped while the sector J, remains 
gapless (which corresponds to the m = 5 plateau phase 
in Sec. lIVp so that the usual spin mode-charge mode sep- 
aration present when m = is no longer valid. 

In the strong-coupling limit of the Hubbard model 
U ^> t and for small hole doping, Hamiltonian (Q]) re- 
duces to the t-J model: 

n = -h E + ^ 

+J\\ E] fthP ' Si+i,p — -£n>i ;P ni+i,p] (5) 

i,p=l,2 

i i,P=l,2 

in which V is the Gutzwiller projector, m^ v = ni iPiCT 
and the same convention is used to label the antiferro- 
magnetic couplings J\\ t ± = 4i| j_/U. In what follows, 
these two microscopic models will be studied numeri- 
cally and the isotropic model will assume t = t\\ = t± 
and ,/ = J|| = Jj_. 

Bosonization allows to study the low-energy proper- 
ties and correlation functions of ID-like systems using 
field theoretical methods. The bosonization Hamilto- 
nians used to describe ladder a 30 i 31 ' 32 i 33 ' 34 fall into two 
classes, coupled bands models to be reviewed in Sec. HVl 
and coupled chains models to be reviewed in Sec.fV] The 
first approach is more appropriate to study the system 
with isotropic parameters. There is however the possi- 
bility, by strongly reducing the interchain hopping am- 
plitude t±, to obtain a cross-over to a regime where the 
coupled chain model is best suited to describe the behav- 
ior of the system (Sec. [Vj. 

III. NON INTERACTING AND 
STRONG-COUPLING LIMITS 

Non interacting system. It is interesting to discuss first 
the non-interacting system using the relations ^ and 



(a) (b) (c) 




FIG. 2: Bands evolution corresponding to FIG. |3] 




FIG. 3: Magnetization curves for the free system with tu = 0.2 
and t± — 1.0 showing a plateau at m = 8 of width 2(t± — 2tii). 



keeping 6 fixed instead of the chemical potential. Two 
main cases are possible: either a strong interband cou- 
pling with t± > 2t\\ or a small interband coupling for 
t± < 2<||. 

In the first case, the only bands which are partially 
filled at low magnetic field are the up and down spins 
bonding bands (0,f) and (0,1) (see FIG. [3 (a)). A first 
critical field corresponds to the complete filling of (0, t) 
(FIG. [2] (b)). This induces a plateau with m = 5 in the 
magnetization curve. This plateau is doping-dependent 
and similar to what has been predicted for other doped 
systemsi 12 i 13 The width of the plateau can be deduced 
from energetic considerations^ 3 - to be 2(t± — 2t») in our 
case. When the plateau ends (FIG. (c)), the band 
(it, t) starts to be partially filled. When the band (0, J.) 
gets empty, all electrons are polarized and the magneti- 
zation curve is constant with value m = 1 — 5. For open 
boundary conditions (OBC), a m = 5 plateau is also 
found as displayed in FIG. [3] (see Sec. HVETI for compu- 
tational details). The width is in good agreement with 
2(t± — 2t»). Note that this plateau does not originate 
from interactions contrary to what will be discussed in 
the strong-coupling approach below and in Sec. IIVE 11 
but is simply due to band-filling effects. 

In the second case, all four bands are partially filled 
at low magnetic field. No plateaus are found but fill- 
ing or emptying successively the bands will induce cusps 
in the magnetization curve and a decrease of the slope 
of the magnetization because fewer electrons contribute 
to the magnetization. First, the band (tt, J.) gets empty, 
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FIG. 4: (Color online) Plateaus in the strong-coupling limit. 
We chose J± = 2.5, Jy = 0.3 and t\\ = t± = 1.0. The m = 5 
plateau is an effect of interactions and the ground-state in 
this plateau has unpaired holes. Still, holes are paired for 
magnetizations with m < m c < S. 




FIG. 5: (Color online) Plateaus phase diagram of a 2-leg lad- 
der for t± — J± — 1.0 and tu = Jy = 0.2. Doping immediately 
destroys the m = plateau present at half-filling while a large 
doping-dependent m = S plateau appears for a small critical 
field. 



next (0, |) becomes completely filled and lastly, (0, |) gets 
empty corresponding to the saturation of the magnetiza- 
tion. An example of such a curve for iy = tj_ with OBC 
is given on FIG. 1161 Two cusps are visible in this curve 
as well as a decrease of the average slope of the magne- 
tization. 

Qualitative remarks on the strong- coupling limit. 
Adding interactions in the system offers the possibility 
to study magnetization plateaus due to interactions and 
also pairing, which does not necessarily mean supercon- 
ductivity. In this limit we have Jy , tu , ij_ <C Jj_ and the 
pairing energy defined as in Eq. (|33| is estimated to be^ 
J± — 2tj_ — 2tu and remains approximately constant un- 
der magnetic field because the average magnetization in 
Eq. is zero. Two ground states are possible at finite 
magnetization: either holes are paired up and all the po- 
larization is due to triplets, or holes pairs are split apart 
and the partner electron on the rung holds the polar- 
ization. The difference between the energies per rung of 




FIG. 6: (Color online) Local densities of holes and spins in 
a system with spins in the ground state of the m = S phase 
computed on a ladder with L = 48 and 6 holes, t± = tu = 1.0 
and J± = 0.5. As J± increases from to Jy , hole pairs form 
in the isotropic limit. Pairing in and above m — S is not 
expected in the strong-coupling limit. 



these two states is: 



m < S 
m > S 



Ae 
Ae 



(Jj_ - 2i ± - 2t||)<5 - Jj_m 
(2Jl - 2t± - 2tu)6-2J 1 _m 



If J±/(t± + in) > 2, a transition is then possible from 
the state with paired holes to the state with unpaired 
holes. The corresponding critical magnetization m c = 
[1 — 2(tj_ + tu)/J±]S is always smaller than 8. This also 
gives a possible scenario for a wide m — 5 plateau due 
to interactions. For large Jj_ and m = 0, we expect the 
ground state to have pairs of holes and singlets mostly on 
rungs so that the system has a large spin gap, which gives 
a m — plateau. However, the m = plateau is much 
smaller than the half-filling spin gap because polarized 
spins in fact gain kinetic energy if they are localized next 
to holes so that it is easier to create them. As the mag- 
netic field is increased, hole pairs are split as suggested 
by the above argument and the partner electron of the 
hole on the rung becomes polarized. When all of them 
are fully polarized, we enter the m = 5 plateau because 
the next spin excitation one can do is to flip a singlet 
on a rung, which cost is approximately J±, The m = 5 
plateau can thus be wider than the m = one. This 
scenario corresponds to FIG.[H where no pairing is found 
in the plateau from the computation of local densities 
of spins and holes with DMRG. Such a scenario could 
be relevant for compounds with lightly coupled dimers 
which could be doped with holes. 

When one starts with unpaired holes at zero magnetic 
field, for instance with tu = J\\ <§C t± = J±, the large half- 
filled spin gap of order Jj_ — Jy is immediately destroyed 
by doping (see FIG. [5]). Increasing the magnetic field fur- 
ther brings quickly the system into a plateau phase with 
unpaired holes. While the magnetization curve is simi- 
lar to the non-interacting one below the m = S plateau, 
its width is clearly controlled by J± — Jy rather than by 
2(i_L-2ti|). 
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An interesting feature on the pairing of holes in the 
m = 5 plateau is the occurrence of a transition when 
J|| is increased from to J± from a state with unpaired 
holes to paired holes (see FIG. [5]) . Having pairing in this 
plateau is a situation that is not expected in the strong- 
coupling limit as we have seen, but for intermediate Jj_, 
we expect that hopping between hole pairs and magnons 
can stabilize hole pairs 2 ^ (see Sec. IIV|) . FIG. [HI suggests 
that a Jii comparable to J± is also needed to have pairing 
in this phase. 



IV. COUPLED BANDS REGIME 

A. Bosonized Hamiltonian 

When the interactions are not too strong compared to 
the interchain hopping, it is reasonable to begin to solve 
the non-interacting band structure and then add interac- 
tions. This is the approach followed in Refs. 32,33 ; 34]. 
The non-interacting band structure is simply formed of 
the bonding band of energy £o,<r(^) ancl an antibonding 
band of energy s n . a (k). The annihilation operators of the 
fermions in these bands are respectively given by: 

1 



1p0,a = -j=(ll>l,a + 1p2,a) 



V2 



(V>l,<r - 1p2,a) 



with 1 and 2 the labels of the chains. In the continuum 
limit, the fermions and 7r are bosonized in terms of the 
boson fields <pk .a (k y = 0, ir) so that we write the fermion 
operator: 



(6) 



in which x = na (a is the lattice spacing) and 



VB,/L,k v ,o 

V27ra 



where R, L are the labels for the right and left movers (see 
FIG. [7]), a is a cutoff (typically a), e R / L = =f1 and the r] 
are Klein factors needed to make the annihilation oper- 
ators of the different fermions species anticommute. We 

+4>R,ky,a] while 

R,k y ,a] = 7T J Hk y ,(T- 

We then introduce the useful bosons <p Ut k with v — c, s 
for charge (resp. spin) corresponding to the symmetric 
(resp. antisymmetric) combination of <pk ,| and <f>k 
The same transformation is performed on the dual fields, 
leading to the Hamiltonian: 



also have the definitions <fik y .<j = \ [4>L,k 
dual fields are 6» fc „,<j = \[4>L,k. 



H 



E 

ft„=0,7. 



dx 
2^ 



[U ky (71-^,^ ) 2 + (d X ^ V . ky ^ 



where [(j>v,k v {x)i'Rv',k' y (rf)] = i8u,v'$k v ,k' y 8( x ~ x ')- In 
the following, we will make the usual approximation 36 
of neglecting the difference between the velocities u k of 
the and 7r bands. This allows us to introduce the linear 
combinations: 

4>c,± = ~^Wc,0 ± 0c,iv), 
4>s,± = -^=(0 S ,O ± ^s.tt), 

with similar definitions for the dual fields Q v ,± = 
7r fU Vt ±, In the most general case, we have to use a 
Z matrix to describe the evolution of the system under 
magnetic field^ When comparing the results with the 
chain models of Sec. [V] it is useful to note that while 
<p c+ = 4> p+ and <p s+ = 4> a+ (we use Greek letters for the 
fields defined in the chain models and Latin letters for 
the fields defined in the band models), while there is no 



simple relation between 



and 



and between 



and 4, 
term: 



The magnetic field couples to the system by a 



H = — / dxd x (j) s+ 

7T 



(7) 



Once interactions are turned on, two types of terms 
appear in the Hamiltonian. The terms of the first 
type are forward scattering interaction terms that are 
quadratic in the fields 4> v .± ■ The terms of the second type 
are backscattering interaction terms, the expressions of 
which were derived in Ref. [33]. Since terms contain- 
ing cos 2(j) s+ cannot appear in the Hamiltonian when the 
magnetization in nonzero, the expression of the backscat- 
tering terms in a magnetized ladder reads: 



ack. 



dx cos 2(9 r _ 



(27ra)- 



cos20 ;j _ 



2g B 
(27ra) 2 



cos 26 s - 
(8) 



From the Hamiltonian (JSJ) , one sees that in the ground 
state the field 9 C - is pinned to (8 C -) = 0. By using the 
results of Ref. [H [Sec. 3.1 and Eqs. (20) and (56)], 
one can argue that in the presence of repulsive interac- 
tions, one must have i4T s _ < 1. Thus, one obtains a 
freezing of the field related to the most relevant oper- 
ator {4> s - ) = f ■ Briefly, when m = S = 0, all fields 
are massive and we have the famous spin-liquid phase, 
often denoted by C0S0. When m = but S ^ 0, the 
system with repulsive interactions is in a C1S0 phase, or 
Luther-Emery phase (LE), with sectors c— and s± being 
massive while the sector c+ is massless, corresponding to 
the charge mode. When m ^ and 5^0, the sector s+ 
becomes massless giving rise to a C1S1 phase. Luttingcr 
exponents associated with these sectors will be denoted 
by K c / S+ . Till now, we have discussed the case of generic 
nonzero magnetization. For the specific case m = 5, how- 
ever, a magnetization plateau can be expected. Following 
Ref. [2a|, we introduce the fields 



J c,± 



V<t>a,±) 



(9) 
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which appear in the term that induces the opening of the 
m = 5 plateaus 



(0,D(0,i) (Jt,D(Jt,l) 



da; cos 2(k FQ + kj 



(10) 



with, from ^ and ([!]), the relation k 



F,0 



7r(n 



am). This term leads to the opening of plateaus when 
n±m € Z as expected in YOA theorem. This condition is 
a commensurability condition which combines both spin 
and charge degrees of freedom. The m = 5 plateau cor- 
responds to the locking of the (f>^ mode. The origin of 
the term (|10l) giving rise to the plateau is the Hubbard 
interaction Un^riii which contains the terms 

(4tM + c It c ot)(4i^1 + 4i c 0i) ; ( n ) 
which, once bosonized, yields Umklapp terms such as 



da; cos 2(k F0 + k Fn )x — 2( 



5 o,T 



VF 



(12) 



B. Superconducting order parameters and most 
divergent fluctuations 

In this section, we define the order parameters for su- 
perconductivity at nonzero magnetization (high field), 
derive their bosonized expressions and deduce their long 
range correlations. Since the SU(2) symmetry is bro- 
ken by the magnetic field, we have to compute the su- 
perconducting correlation functions (A^, (x)A^, (0)) in 
various channels A. We use the following microscopic 
definitions for the pairing operators A^, (n) at rung n: 



singlet: A|^ 



(n) = 2J vc rr7 c r > 



triplet: \ AL (n) = c r jC r 'j 



(13) 
(14) 



with r — (n, 1) and r' = (n, 2) for next-nearest neighbor 
pairs created on a rung and r — (n,p) and r' = (n+l,p) if 
created on the leg p. Contrarily to the case of the coupled 
chain regime (Sec. [Vj where the Fermi wave-vectors are 
the same in both chains, in the case of the band regime 
the Fermi wave-vectors of the two bands are different 
(k FQ 7^ k F7T ) as can be seen on FIG. [I] As a result, 
a more detailed derivation of the bosonized expressions 
starting from lattice expressions in a two-chain Hubbard 
model becomes necessary. 



1. Bosonized form 

Using the bosonized form ([6]) of the fermion operators, 
we can express these order parameters as a function of 
products of the ipR/L,k I n this section, only compo- 
nents with the dominant contribution will be kept, i.e. we 




e(k)A 



(71,1) (7t,t) (0,1) (0,1) 




FIG. 7: Schematic representation of bosonized bands in the 
coupled band models. 



will neglect terms of the form ^R^R an( i tyitPL- These 
dominant contributions correspond to pairing with the 
lowest total momentum q. The case of 2k F triplet pairing 
will be discussed in Sec. IIV B 21 At finite magnetization, 
the equality of the velocities and equation Q ensures 
that the lowest momentum is q — k FQ , — k F0 ^ = nm. 
To simplify the expression of the operators we will use 
extensively the following results on a pinned field ip (for 
which (ip) = cste): we can replace (/(</?)) by /((</?) ) 
and the dual field has exponentially decaying correlation 
functions 

Starting with the interband order parameters, which 
read: 



i»Rfi,atyh,-K,o 



--0c- -f 



'.+ -«.-)] 

) 

.+ -*-)] 



we note that they are all proportional to e 1 ^"- and thus, 
since 9 C - is pinned, their correlations decay exponen- 
tially. In other words, power law decay is possible only 
for intraband superconducting correlations. 

The intraband superconducting order parameters in 
the three spin channels read respectively: 
-Intraband singlet: 



+ +0 c _- CT (0 s+ +0 3 „)] 
+-*.-)] 



— 6f c _— cr 



y^^^R,0,a^L,Q~a ~ ° 
<t a 

cr cr 

-Intraband triplet S z = 0: 

cr cr 



-Intraband triplet S z = 1: 

1pR,0,tT^L,0,tT ~ < 
1pR,ir,alf)L,ir,cr ~ « 



i[9 c+ +9 c -+<T(e s+ +e,-)] 
i[6i 0+ -ec-+o-(9a + -e»_)] 



To determine which forms of superconductivity will be 
dominant, we need to express the leg and rung order 
parameters in terms of the intraband order parameters. 
We have for the rung singlet order parameter, 



IpR^^L^-a), 

(15) 
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Whereas the rung triplet S z = is given by: 

A U( X ) ~ E ^e lCT<?a; (V'fl^ :( T'0L : O -a ~ 1pR,0,-<rll)L,K,o), 

therefore, since it is composed only of interband terms, 
we expect that its correlation will present exponential de- 
cay. Lastly, for the rung triplet S z = 1 order parameter, 
we find: 

A TT( X ) ~ ^(^,0,1^,0,1 - •tpR,TT,li>L,Tr,t), (17) 

Turning to the leg singlet order parameter, we find that 
it reads: 



'■° a ipR,o,aipLfl -<r + (0 -> tt) 



(18) 

Note that because we can neglect interband coupling, the 
expression is the same on both legs. Thus leg-leg correla- 
tions on the same leg and between the legs will have the 
same sign as found numerically in FIG. Analogously, 
the leg triplet S z = operator reads: 

A U X ) ~ E [MkF, a)i>R,o,vi>L, ,-* + (0 -► tt)] 

(19) 

Note that if we take the limit a — ► this term disappears 
as happened for the rung triplet. Finally, the leg triplet 
S z = 1 reads: 



AUx) 



'■° a 4>RAAipL,o,i + (0 -> tt) 



(20) 



Since < 1, we have (</>«-) = f ■ Then, all the in- 
traband triplet S z — 1 have exponentially decaying cor- 
relations. As a consequence, both the leg and the rung 
triplets with S z = 1 have exponentially decaying corre- 
lations. The behavior of these correlations with higher 
2fcp momentum will be discussed in Sec. lIVB~2l Another 
consequence of the ordering of the field (j> s - is that: 



i>R,0,<r1pL,Q,- 



gif <? e i(8c+-cr<t> ! ,+ ) ^ 



Provided that Uf,o ^ kF,*, the leg triplet S z = or- 
der parameter does not vanish. We find that it has the 
same critical exponent as the rung singlet order parame- 
ter, namely ^(K^ + K s+ ). It is larger than the m = 
critical exponent 1/(2K C+ ) because of the appearance of 
the massless boson 6 S+ . The fact that we have a finite 
momentum pairing with q = nm which is very likely to 
occur in a one-dimcnsional-like system and is the signa- 
ture of the FFLO mechanism. 



the band representation, it is indeed possible to find a 
component of the rung-rung triplet correlations that has 
power law decay. Let us consider the next terms in the 
expansion of the rung triplet operator: 

In bosonized form, this operator reads: 

A^ (j2fc (x) ~ e i ( k F,0+ k F,^) X e i l e o+-<Pa++Cr{0 s+ -<p s+ )] 

It has power law correlations with a critical exponent 
1 



K c+ + K c l + A' s+ + Kj 



and, from equations ((3]) and ([4]), the associated wave- 
vector is k F0 + k Fn = ix(n + am). Note that its critical 
exponent is always larger than the one of the rung singlet 
order parameter. 



3. Charge- density waves order parameters 

Here, we address the question of the CDW exponents 
under magnetic field. First, the 2kp CDW order param- 
eter rii, which vanishes exponentially at zero magnetiza- 
tion, contains terms such as: 



f> c -+<7(t 



-r(e c _+<re o _)] 



with r = ± for R, L and a wave- vector k F7r — k F0 . These 
order parameters decay exponentially because they are 
proportional to e l ^ c - . A last possible term is 



+ -e c _+<T(</, s+ -0 s _)] 



with a "2/cp" wave vector 7r(n + am). It decays expo- 
nentially because of e s ~ . Finally, no 2fci? CDW are 
expected under magnetic field. 

Secondly, the 4fcp CDW order parameter— is rof and 
has an exponent 2K C + at zero magnetization. This order 
parameter contains terms such as 

which have a "4fci?" wave vector 2ix(n + m) and, for 
nonzero magnetization, an exponent 2{K C+ + K s+ ). We 
also have terms like 



2. S z = 1 Triplet with a 2/cf momentum 

Another notable result is the presence of rung-rung 
triplet S z — 1 correlations with 2fcp oscillations. Within 



with a wave-vector 2{k F + kp^) and an exponent 2K C+ 
which is not affected by the magnetic field. This last term 
can compete with the superconducting order parameter 
to be the most diverging fluctuations depending on K c+ 
and K s+ . 
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Type of operator 


fermion expression 


dimension 


wave- vector 


S z = 1 triplet 


1pR0li>L0L 
4>RKl4>L07 


4K7 ~*K2 2K3 






Fir ™F0 


S z = triplet or singlet 


IpROflpLOl 
lpRTr1*pL0l 


^(^r 1 + ^2) + ^(if 2 - 1 + ^1) + 1(^3 + k^ 1 ) 

^(^r 1 + + '^(K- 1 + K{) + UKs + K- 1 ) 


k^ 

^FO ^FO 
k^ — k^ 



TABLE I: The superconducting operators for the Hubbard model with an empty band corresponding to the coupled band 
regime just above H c . 



Type of operator 


fermion expression 


dimension 


wave-vector 


SDW'/CDW 


4> RO l1pL01 


Ki , Ki , K a 
4 ' 4 ~ 2 

Ki , if 2 , K 3 
4 ' 4 "T" 2 

2^2 
^1 1 Ti 1 1 
4 ' 4 ' 2if., 


^ K FtV 

-4- i- T 


SBW X ' V 


^RO^LOL 


^(ifr 1 + K 2 ) + 2±^(^ 2 " 1 + K X ) + ±(K 3 + K3- 1 ) 
^(ifr 1 + AT 2 ) + 2±^I(^ 2 - 1 + Ki) + §(JT 3 + K- 1 ) 


i-T -1- k^ 
"•FO ' ™F0 

k^ -4- k^ 
"■Fir ">F0 



TABLE II: The CDW/SDW operators for the Hubbard model with an empty band corresponding to the coupled band regime 
just above H c . 



Interpreting the superconducting critical Held 
H c as a band-filling transition 



In this section, we identify the observed superconduct- 
ing upper critical field H c with a band-filling transition. 
Such a transition will ungap all 3 remaining sectors, lead- 
ing to an enhancement of the correlations exponents. We 
describe the system just above the transition and com- 
pute the various possible exponents using bosonization. 
We consider a case where the Fermi energy is positioned 
in such a way that the band (ir, J.) is empty while the 
three other bands remain partially filled (see FIG. Q] for 
illustration). The band (n, J.) being empty has important 
consequences. Projecting out the high energy subspace 
where the band (0, J.) is occupied by a single electron, 
one gets in lowest order: 



t 1 f 



(21) 



and as a result the on-site Hubbard interaction reduces 
to: 



(22) 



completely gapless. The bosonized Hamiltonian reads: 

(23) 



H 

Ho 

Hij 



Ho + Hi 



«e{(o,t),((U).(f.T)} 
Ua 
2^ 



dx 



2 f dxd x (f> 0i i(d x (j)o : -\ + d x cj) v ^) 



(24) 



This Hamiltonian can be fully diagonalized. In the 
general case, where the Fermi velocities are all different, 
one needs first to perform a rescaling: H, — > ^Ju/vf,v^-v 
and (f> u — > § v j \JujvF. v where u is an arbitrary quantity 
with the dimension of a velocity and then diagonalize the 
matrix: 



-'.f.ot 



Ua 







^y/V F ,0]V F fil 
V F,0l 



\ 



J$y/VF,KtVF,Ol 



(25) 

Afterwards, one performs the inverse rescaling on the 
diagonal matrix to obtain the velocities of the modes and 
the associated Luttinger exponents. To simplify the al- 
gebra, we will assume that all the Fermi velocities are 
equal. Then, we find that the following combination of 
fields 



In the absence of commensuration between the differ- 
ent bands, the resulting bosonized Hamiltonian has only 
forward scattering interactions. Thus, its excitations are 







( 00,1 \ 


) 






V / 



(26) 
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FIG. 8: (Color online) Absolute value of superconducting cor- 
relations for zero and a finite magnetization in an isotropic 
doped ladder with 5 = 0.063. Notation are "Rg" for rung, 
"Lg" for leg, "~ 0" for numerically irrelevant signal, "exp" 
for exponentially decaying correlations and "7rm" and "2fcF" 
for the wave- vectors. Note that the S" — 1 "2fcF" oscillations 
are smoothed out by taking the absolute value. M = 1000 
states were kept with discarded weight of order 10 -6 . 



diagonalizes the interaction, and: 



U\K\ = u 2 K 2 = u 3 K 3 = v F 
Ui Ua 

u 2 
K 2 

K 3 



VF 
V F 
Vp 



2tt 2 V2 
Ua 



(27) 
(28) 

(29) 
(30) 



In terms of these new fields, the fermion operators read: 



VrfiJ i[i(ei-r<t> 1 )-^(e2-r<f> 2 )+^ s (e 3 -r4> 3 j\ 

\j2na 

, _ r)r,n,T ^[i(8 1 - r 1 )-i(fl 2 -^ 2 )--^(e 3 - r 3 )] 



From these expressions, it is possible to obtain the vari- 
ous superconducting order parameters and their critical 
exponents. The results are gathered in table Q] and ta- 
bic |TT] on which one must note that the correlation expo- 
nent will be twice the dimension of the operator. 



D. Numerical results on superconducting 
correlation functions 



We have computed superconducting correlation func- 
tions with DMRG using the definitions of Eqs. (fl~3"|) and 




32 r 48 



FIG. 9: (Color online) Same correlations as FIG.|8]but using 
the Hubbard model with 6 = 0.063 and U/t = 8.0. M = 1600 
states were kept with discarded weight of order 10 -6 . 



(|T4|) . As predicted by bosonization results, only modified 
d-wave singlet superconductivity is found when m = 
(rung-leg correlations have an opposite sign to rung-rung 
and leg-leg correlations, see Ref. 16] for a bosonization 
discussion), while S z = leg-leg triplet and S z = 1 rung- 
rung triplet correlations emerge under high magnetic field 
(see FIG. U for the t-J model and FIG. M for the repul- 
sive Hubbard model). Furthermore, these correlations 
oscillate with wave- vectors q = nm for the S z = chan- 
nel and q — 2kp for the S z = 1 channel. The q = nm 
relation has been checked in all the FFLO phase (see 
FIGfTTl and discussion in Sec. IIVE]) . by fitting the oscil- 
lating correlations with Acos(qr + (p)/r a (see FIG HOp . 
confirming the FFLO-like mechanism. We observe that 
S z = rung-rung triplet correlations have an exponential 
behavior while rung-rung singlet are algebraic, in agree- 
ment with the prediction from bosonization. Since from 
the point of view of the rotation symmetry these two or- 
der parameter transform in the same way, the difference 
must be caused by their different transformation prop- 
erties under interchange of the legs. Note that because 
of the sm(k F k a) factors in the leg singlet and S z = 
leg triplet, the cancellation of the intraband terms ob- 
tained in the case of rung order parameters is absent. 
This explains the observation of power-law correlations. 
Concerning the orbital part of pairs, singlet Cooper pairs 
have a mixed s-wave and (i-wave structure and S z = 
triplet can be considered to first approximation as p-wave 
pairs along the legs (with a symmetric superposition of 
the two legs). 

The good agreement with bosonization predictions re- 
lies on the fact that in the doped isotropic t-J ladder, 4 
Fermi points exist with approximatively 35 fci?o ~ 37r/5 
and kF,Tr ~ 27r/5 when m = and that the assump- 
tion of equal Fermi velocities is numerically reasonable at 
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o singlet 

□ S z =0 triplet 
— q = ran 




from L = 64 
5 = 0.0625 
J/t = 0.5 



0.15 

m 

FIG. 10: (Color online) Verification of the q = nm relation 
of FFLO-like pairing at finite magnetization in an isotropic 
ladder with J/t = 0.5. 



low hole doping. The existence of a small leg-leg triplet 
component (about 100 times smaller than the rung rung 
triplet) in the two leg t-J ladder considered may therefore 
be explained by the persistence of a small band splitting 
in this strongly coupled ladder system. 

We now would like to compare precisely the algebraic 
decay exponents a s and a* of the singlet and S z = 
leg-leg triplet correlations. DMRG is known to often un- 
derestimate correlation functions 25 for a fixed number of 
state kept M. In order to capture the behavior in the 
thermodynamic limit, we computed the correlation func- 
tions for fixed M ranging from 800 to 2000 and lengths 
L = 32, 64, 96 and 128 (we worked at fixed 5 = 1/16 
and m = 3/32 so that these are the only accessible clus- 
ter sizes). We extracted a(M,L) by fitting the data. 
Then, we can extrapolate a in the M — > oo limit to get 
a correct exponent a(L) at size L from which we can 
do a finite size scaling. On FIG. QT] we see that for a 
given size, a(M) decreases as M is enlarged, roughly like 
a(M) ~ 1/M (see FIG. [13] (b)) as was previously found. 
We note from FIG. [T2] that the larger the system, the 
larger M is needed to reach a good convergence. Con- 
vergence as a function of the discarded weight is also 
given on FIG. [13] and has qualitatively the same behav- 
ior. For L = 128, the convergence is slower with M than 
for L = 96, probably because we would need larger M 
to have a correct accuracy (the 1/M might be realized 
for large enough M) so that we believe the results are 
not as reliable as for L — 64 and 96 (larger M would 
be very expensive numerically). On one side of FIG I13I 
(c), M = 2000 is too small (for L = 128) while on the 
other side, it is difficult to extract a because the system 
is too small to resolve enough oscillations (when L = 32). 
Extrapolations can be tentatively done with some uncer- 
tainties which can be roughly estimated. Therefore, we 
can infer from FIGH3] that a s = 1.54 ± 0.15 is greater 
than a 1 — 1.17 ± 0.15 which pleads in favor of dominat- 



M = 400 

600 

800 

1000 

1200 

1400 

1600 




FIG. 11: (Color online) Convergence with the number of kept 
states M for a fixed length L = 64. The plotted super- 
conducting correlation is the singlet one and parameters are 
m = 0.094, 8 — 0.063 and J/t = 0.5 which corresponds to the 
phase above the m = 5 plateau and below the superconduct- 
ing upper critical field H c . 
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FIG. 12: (Color online) Convergence with size L for the same 
parameters as in FIG. 1111 M being fixed. 



ing S z = triplet correlations in this part of the phase 
diagram (above the m = 5 plateau and below H c ) but the 
difference is rather small. Finite size effects could explain 
that the bosonization prediction a* = a s is not realized 
on numerical results. Note that a < 2 is sufficient to have 
a divergent superconducting susceptibility in the system. 
Concerning CDW exponents, they are difficult to com- 
pute numerically because of OBC but we will see that we 
can deduce from the two-particle gap of FIG. [TS]that the 
system is in a superconducting phase. 
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FIG. 13: (Color online) Extrapolation of fitted exponent as a 
function of discarded weight w (a) and of the inverse of the 
number of kept states M (b) for different lengths L. Extrap- 
olated results from (b) are tentatively extrapolated vs L in 
(c). Large error bars occur when L is too small (for L = 32) 
and because M is too small (for L — 128). Parameters are 
given in the caption of FIG I11I 
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FIG. 14: (Color online) Normalized 4-points spin correlations 
T(x) = (S2 {x)S} (x)S2 (0)Si (0)} for various magnetization 
in the isotropic t-J ladder with J/t = 0.75 and S = 0.063. 
Their decay exponent is 2K~1 which gives access to K a +(H). 
Correlations were computed on a system with L = 64 and 
M = 1200 states kept. 



Evolution of the Luttinger parameter K a + 

We now want to have access to the evolution of the 
Luttinger exponent K s+ of the gapless mode which ap- 
pears at finite magnetization. We remark that the 
bosonized form of the triplet creation operator on a rung 
S2 {x)Sf (x) contains dominant terms such as 



-i2[0 s _+e s+ ] 



(31) 




•-•J/t 


= 0.35 


■-■J/t 


= 0.5 


♦-♦J/t 


= 0.75 
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1.5 2 

H/A 



2.5 



FIG. 15: (Color online) K a+ (H/A a ) from DMRG computa- 
tions on a system with parameters of FIG. 1141 



giving an exponent 2K 8 } since {(f> s -) = tt/2. The wave- 
vector associated with this operator is (kp Q 



(k] 



'4J 



since, if the difference between the Fermi 



velocities of the and ir bands is negligible, we have 
jm. Note that 2if s+ 1 is the smaller ex- 



it 1 

F,k y 



k l 

K F,k y 



ponent for terms with this wave-vector. Because K s+ < 
1, we could expect smaller possible exponents such as 



2K» 



or 



K 



For the first, one can show 
or 2\<t> s+ + ctO c -] cannot ap- 



that fields 2[<f> s+ ± 

pear in the decomposition of S^S^. For the other, fields 
cj> s + ± 9 S + ± 9 C - ± 4> s - cannot be decomposed in terms of 
the right and left movers fields with factors ±. Finally, 
the decay exponent of the correlations of this order pa- 
rameter is expected to be 2 and depends only on K s+ 
and not on K c+ . Numerically, we have computed the 
correlation function (5'^"(a;)5'J , "(a:)5'^"(0)S' 1 "(0)). Data for 
the isotropic t-J model are displayed on FIG.fLlland show 
a power-law decay with an exponent larger than 2 and 
a wave- vector approaching zero. Points are given as a 
function of H/A s , A s being the finite-size spin gap of the 
system. We thus have access to the evolution of K s+ as a 
function of the magnetic field (see FIG. IT5|) . The general 
behavior is that K s+ < 1 and increases with magnetic 
field towards the limit K s+ — 1 at high fields. Similarly, 
the Luttinger parameter K s+ has been characterized us- 
ing Bethe ansatz in the context of the SO (8) descrip- 
tion of two-leg half-filled Hubbard ladders-- which gave 
K s+ < 1. When doped, the SO (6) description^ gives 
the same constraint (results not shown) as well as an in- 
crease of K s+ at high magnetic field. However, a detailed 
comparison of SO (6) predictions and numerics would re- 
quire to determine first the conditions and parameters 
at which both approaches match and this is beyond the 
scope of this article. Lastly, the values of K s+ obtained 
numerically should be slightly larger in the thermody- 
namic limit than what is found because of finite size and 
finite M effects as explained above. 
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FIG. 16: (Color online) Magnetization curves for the isotropic 
t-J, Hubbard and non-interacting systems. Irrational plateaus 
are clearly visible at m = S. H c represents the superconduct- 
ing critical field. The energy scale J is defined as 4t 2 /U in the 
Hubbard model. The non-interacting system curve has been 
rescaled (H — > H/4) for clarity. 
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5 

FIG. 17: Phase diagram for the isotropic two-leg t-J ladder 
in the (H, S) plane for J/t = 0.35. The large FFLO phase is 
delimited by the upper limit of the m — plateau and the 
critical field H c . Note also that in this part of the diagram, 
the m = 5 phase is not expected to be a superconducting 
phase but rather metallic. 



E. Generic phase diagram of the t-J model 

In this section, we discuss the generic phase diagram in 
the (H, S) plane for the isotropic t-J model on the basis 
of DMRG and bosonization results. 



1. Doping- dependent magnetization plateaus 

The magnetization curves m(H) of FIG. HH obtained 
within the Hubbard and t-J models, display plateaus for 
to = and m = 5. Energies E(rih,S z ) were computed 
keeping M — 1600 states with the single-site method 
proposed by White^ 2 - (we used a noise level of 10~ 6 ) at 
fixed hole number rih and total magnetization S z . Mag- 
netic fields are deduced using H(S' Z ) = E(rih, S z + 1) — 
E(n h ,S z ), and interpolated with \H{S X ) + U(S Z - 1)]/2 if 
they do not belong to a plateau. The m = plateau sim- 
ply corresponds to the well-known spin-gap of the doped 
ladder. The m = S plateau, exists at small doping for 
continuous values of 5 (see FIG.[l"7land Ref. (23|) and thus 
falls into the classes of doping-dependent irrational mag- 
netization plateaus predicte d 12 ' 26 by Cabra et at. It can 
be understood as a Commensurate-Incommensurate^ 
(C-IC) transition so that we expect the critical exponent 
of the magnetization as a function of the magnetic field 
to be 1/2. In the plateau phase, the mode cjx^ is locked 
but the sector (J,, +) remains gapless, leading to a metal- 
lic phase. Note that (j>t is a superposition of both spin 
and charge modes so that charge and spin are no more 
independent modes in the plateau phase. 

The Hubbard and t-J models give qualitatively the 
same behavior at low magnetization. The spin gap 



(to = plateau) in the Hubbard model with U/t = 8 is 
about half of the spin gap of the t-J model with J/t = 0.5 
as it was previously found^ but the irrational plateau 
has roughly the same width. A larger U gives a slightly 
larger plateau. Both models displays a singularity of the 
magnetization curve, more or less pronounced, near the 
superconducting critical held H c . The location of H c 
in FIG. [THl have been roughly determined by looking at 
local densities (see Ref. |23j for method) since Cooper 
pairs break down above this field. For larger U, the 
superconducting critical field is smaller. These results 
are also observed in the t-J model having in mind that 
J/t ~ t/U (see Sec.|PVE2Jand FIG. US). For higher mag- 
netization states, Gutzwiller projection induces a rather 
different behavior between strongly interacting systems 
(Hubbard U = 16 and t-J) the Hubbard model with 
U = 8 which displays the expected square-root-like be- 
haviors near critical fields in good agreement with the 
band-emptying scenario proposed in Sec . 11 V 01 At last, a 
quick comparison with the isotropic non-interacting sys- 
tem proves the non-trivial role of interactions in the ap- 
parition of these plateaus. 

The magnetic part of the phase diagram of the t-J 
model with J/t = 0.35 has been computed (see FIG. ITTj) 
and is very similar to the one obtained with J/t = 0.5 
(see Ref. [2j|) but have larger to = S plateaus. We pro- 
pose that such a phase diagram is generic for the isotropic 
t-J and Hubbard models under Zccman effect at low dop- 
ing and for parameters 0.25 < (J/t ~ 4</C7) < 1.0 cor- 
responding to the strong-coupling regime. Varying J/t 
modifies the width of the different phases as proposed in 
FIG. 1191 Lastly, we note that 5=1/4 corresponds to the 
end of the plateau both for J/t = 0.35 and J/t = 0.5. 
Furthermore, the magnetization curves are strongly mod- 
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FIG. 18: (Color online) Two-particle gap (a) and pairing 
energy (b) for isotropic ladders as a function of magneti- 
zation. An anomaly in the two-particle gap at m = 8 is 
clearly visible while the pairing energy remains finite up to 
the superconducting-metallic state transition. m c represents 
the magnetization corresponding to the superconducting crit- 
ical field H c - 

ified with 5 > 1/4 (data not shown) which suggests that 
this hole density corresponds to a singular point. 

2. Charge gaps and exceeding of the Pauli limit 

To study the charge degree of freedom of the system, 
we computed the two-particle gap A2 P (related to the 
inverse compressibility) and pairing energy A p as a func- 
tion of the magnetization m to discuss the nature of the 
ground state. With rih the number of holes and S z the 
total spin along the magnetic field, we have the defini- 
tions 

A 2p (S z ) = E(n h + 2,S Z )+E(n h -2,S Z ) 

-2E(n h ,S z ), (32) 

and 

A P (S*) - E(n h -1,S Z + 1/2) + E(n h -1,S Z - 1/2) 
-E(n h ,S z )-E(n h -2,S z ). (33) 

The evolution of these gaps under increasing magneti- 
zation is displayed on FIG I18I Superconducting correla- 
tions are also algebraic when the pairing energy is finite 
from FIG. 6 of Ref. (23[. These observations confirm that 
the system is in a superconducting state, except in the 
m = 5 plateau where a finite two-particle charge gap is 
found as well as a strong reduction of all superconducting 
correlations from FIG. 6 of Ref. [23j]. However, the field 
(f>~^ remains gapless in the irrational plateau while other 
fields are gapful. This suggests that the system is in a 



metallic phase in this plateau phase. Increasing the mag- 
netic field brings the system back into a superconducting 
phase, giving an example of reentrant superconductiv- 
ity which shares similarity with a proposal for a bilayer 
system.— The finite value of A2 P in the thermodynamic 
limit within the plateau phase can be inferred by us- 
ing Eq. (132|) and the constraint m — S for continuous 
values of hole doping, which are more general hypoth- 
esis than the particular case under study. Indeed, it is 
straightforward to find that it exactly equals the width 
E(nh, S z + 1) — E(rih,S z — 1) of the plateau in the ther- 
modynamical limit. This is coherent with the numerical 
results of FIG. [18] obtained at fixed 6. It simply means 
that removing a hole pair in a system which has m = 5 
(i.e. the same number of hole pairs and magnons) is ener- 
getically equivalent to adding a magnon in a system with 
m! = 8' at a slightly lower density. Thus, one has to pay 
the energy gap of the ml — 5' plateau for that, and this 
gap equals the m = 6 gap in the thermodynamical limit. 

We now come to the computation of the theoretical 
Pauli limit H p and of the actual superconducting critical 
field H c . The Pauli limit is defined by equaling the con- 
densation energy at T = and H = to the energy of 
the electrons coming from a broken pair and stabilized 
by Zeeman effect X We have access numerically to the 
condensation energy of a pair of electrons et zero mag- 
netization through A p (S' z = 0). Paired electrons in a 
singlet state do not take advantage of Zeeman effect con- 
trary to unpaired electrons which can be polarized along 
the magnetic field. The total energy of these two elec- 
trons is 2E(nh — 1, +1/2) — 2 x i x H, to be compared 
with the total energy of the same system with paired 
electrons, E(n h , 0) + E(n h - 2,0). By looking at Eq. ([33]) 
with S z — 0, we find that the Pauli field simply reads 
Hp = A p (,5' 2 = 0). What actually occurs in the system is 
a transition to a FFLO state which changes the nature of 
the ground-state so that the Pauli limit can be exceeded. 
The superconducting upper critical field is deduced from 
the location at which A p (S z ) crosses zero which gives a 
critical magnetization m c and in turn the critical field H c 
reported on FIG. [17] It is important to note that these 
calculations are exact and with no approximation. The 
superconducting critical field H c can also be interpreted 
in the bosonization language as the emptying of the (n, J.) 
band (see Sec. HVC]) . From FIG[l7]and FIGEES it is clear 
that Pauli limit is exceeded at low doping and for a wide 
region of parameters. This exceeding can be associated 
with the FFLO mechanism discussed in Sec. IIVDI from 
the behavior of the superconducting correlations. 

At last, we remark from FIG. [18] that the pairing en- 
ergy has a discontinuity upon magnetizing the ladder (in 
the m — ► limit) very similar to the known discontinuity 
of the spin gap 46 i 47 upon doping the ladder (in the 5 — > 
limit). Note also that the strong reduction of the pair- 
ing energy and of the spin gap for 6 — 1/4 observed on 
FIG. [T7] is related to the proximity of the CDW phase 16 
which occurs for J/t ~ 0.2. 
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FIG. 19: (Color online) H c , H p and m = 8 plateau on an 
isotropic two-leg t-J ladder doped with 2 holes as a function 
of J/t (this corresponds to the 8 — > limit of the phase di- 
agram of FIG. 117) 1 . The red curve (squares) corresponds to 
the upper limit of the m = plateau and the lower limit of 
the m = 8 plateau while the black curve (circles) to the upper 
limit of the m = 8 plateau. The FFLO phase is delimited by 
the black curve (circles) and the blue curve (triangles). Re- 
sults are extrapolated from DMRG results on systems with 
L — 32, 48, 64. For small J/t, a Nagaoka phase is found but 
because of strong finite size effects, we do not discuss this part 
of the diagram. Lines are guide for the eyes. 



3. Phase diagram 

The phase diagram of FIG. [I7]is proposed to be generic 
in the coupled band regime for strong repulsive interac- 
tions. In addition to the magnetic properties described 
above, the system is in a LE superconducting state in 
the m = plateau, in a Luttinger Liquid superconduct- 
ing phase below H c and in a metallic phase in the m = 8 
plateau phase. Above H c , the system switches to a metal- 
lic phase with three ungapped sectors which couple both 
spin and charge degrees of freedom. The central charge 
c is expected to take the following values in each phase 
encountered as the magnetic field is increased (taking 
for example a vertical cut on FIG. [17] with 8 — 0.063): 
c = 1 (to = plateau), 2 (FFLO below the m = 8 
plateau), 1 (to = 5 plateau), 2 (FFLO above the m = 8 
plateau), 3 (above H c ), 2 (saturation phase). Note that 
one could expect a last transition c = 2 — * 1 for certain 
parameters, corresponding to a situation where the (0, f ) 
is completely filled. Since such a transition would induce 
a cusp in the magnetization curve and there is no sign of 
it on FIG. [TBI we conclude that c = 2 is the value for the 
saturation in this regime. 

To evaluate the role of J/t on this generic phase di- 
agram, we computed various energy gaps and critical 
fields in the <5 — limit, i.e. for two holes on a lad- 
der with L — > oo. Results are extrapolated from systems 
of length L = 32, 48, 64 and are displayed on FIG. QU 
First, we note that for J/t very small, a Nagaoka phase 



competes the LE phase and induces strong finite size ef- 
fects in the ground states, thus the displayed boundary 
between these phases does not correspond to the ther- 
modynamical transition. Hence, we only focus on points 
for which the LE phase is stable (on finite size systems) 
which are found to be qualitatively for 0.25 < J/t < 1.0. 
For large J/t the model undergoes a phase separation 
between holes and spins and we are not interested in the 
behavior near this instability. In this region of param- 
eters, we have the following behavior: the width of the 
plateau increases for decreasing J/t (for J/t = 1.0, the 
to = 8 plateau is hardly visible at finite density, data not 
shown), while the exceeding of the Pauli limit increases 
with J / 1 (which simply follows the increase of the pairing 
energy). 

These elementary excitation gaps at 8 — > + can be 
related to the dynamics of the system^ 7 - when m = 0. 
The difference between min(A p ,AM) and the spin gap 
corresponds to the binding energy of the resonant mag- 
netic mode. It is interesting to note that the upper limit 
of the m = 8 plateau is approximatively independent of 
J/t and corresponds to a free magnon (which is nothing 
but the spin gap of the undoped system Am ~ <//2). 
Indeed, once the hole pair is bound to a magnon, the 
next magnetic excitation is to create a magnon in the re- 
maining undoped background, which cost is Am, slightly 
renormalized by the scattering with the holes-magnon 
bound-state. This supports the phenomenological mech- 
anism proposed in Ref. [23| for the opening of the ir- 
rational plateau. The binding of the hole pair to the 
magnon comes from the gain in kinetic energy holes have 
in a locally ferromagnetic environment. It is thus ex- 
pected that the binding energy increases with t/J. It is 
also likely to find such a bound state in the vicinity of 
a Nagaoka phase which here competes the LE phase. It 
also explains the decreasing of the pairing energies as J / 1 
is reduced. Note that previous studies of this resonant 
magnetic mode with exact diagonalization 47 and periodic 
boundary conditions agree with these results, and hence 
suggest that open boundary conditions do not affect this 
generic phase diagram. 



V. COUPLED CHAIN REGIME 

A. Dominant exchange model 

By strongly reducing the interchain hopping ampli- 
tude, one can reach the coupled chain regime in which 
interchain hopping is an irrelevant perturbatio n 48 i i 50 i 51 
and the relevant interchain couplings are either the 
Josephson coupling (for attractive intrachain interac- 
tions) or the exchange coupling (for repulsive interchain 
interaction). An important point is that even if the bare 
model has only interchain hopping, Josephson and ex- 
change couplings are generated by the Renormalization 
Group (RG) flo w 49 i 50 i 51 and, as a result, the effective 
model always contain this type of interactions. In the 
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rest of the paper, we will call this strong coupling limit 
the "chain representation" . The bosonizcd Hamiltonian 
describing the two uncoupled chains reads^: 



tion of a magnetic field strong enough to induce a C-IC 
transition^ 



H 



Ef dx 
7 2^ 



(34) 



Dominant exchange 



where we have dropped terms oc gx± cos V8</>cr,i since 
these terms are marginally irrelevant in the case of re- 
pulsive interactions. In the case of attractive inter- 
actions, they also become irrelevant upon the applica- 



For small J\\/tn ratios, the analysis of the scaling di- 
mensions in the t-J model on a single chair*^ shows that 
the dominant interchain coupling is the exchange one. 
This term reads 



J 



2J ± a 
(27ra) ; 



cos 



V2(<f> P ,i - 4> p ,2) 



cos V2(9 aA - ff)2 ) + ^ cos V2(<f><r,i - <t>a.i) + \ cos V2(0ff,i + <f>tr,2 + 2mnx) 



(35) 



where we have the usual definition 9 v< i — J ifO- v ,i- For 
nonzero magnetization, the last term is oscillating and 
drops from the Hamiltonian. It is convenient to introduce 
the new fields 



1 



v,± = — 7=(0v,l ± 0v,2), 



With this transformation, the exchange term is rewritten 
as 



2Jj_ 
(27ra) ; 



dx cos 2(j) p 



cos 29 a - + i cos 2(f> „ 



(36) 



Moreover, the two chains being equivalent, the chain 
Hamiltonian is rewritten as 



r=± 

U=p,P 



u u K v {-kW v ^) 2 + -^-(d x <fi vr ) 2 



(37) 



Obviously, the Hamiltonians describing the fields <fi a + 
and 4> p + are purely quadratic indicating that the total 
charge and the total spin excitations are gapless. The 



fields 



and 



arc described by a generalized sine 



Gordon model. The sine Gordon interaction term 
can be treated within a RG analysis. The scaling di- 
mensions of the term cos 2</> p _ cos 20 CT _ is K p - + K~_ 
and scaling dimensions of the term cos2^_ cos20o-_ is 
K p - + K„- . From the analysis of Ref. [30| we can con- 
clude that in the case of interest we have K a - > 1 and 
cos 2<f) p - cos 26 a - is the most relevant term and for anti- 
ferromagnetic J± the ground state has (<fi p -) = (mod 
7r) and {6 a -) = f (mod n). 



2. Most divergent superconducting fluctuations 

The expression of the intrachain order parameters in 
terms of the fields in ||3"T|) can be found in Ref. [3£| . When 
re-expressed in terms of the ± fields and denoting q — 
Trm, the singlet operator reads: 



O 



SS,i 



(38) 



E 



and, for triplets operators, we have: 



O 



TS,i,S*=l 



giqcx e i{ e P+-(~Y e P-) e - i<7 i 

(7 



(39) 

1 

(40) 



7+ -{-ye a -)_ 

For the interchain operators, the singlet operator reads: 



ss 



E 

a 

E 



(TC 



iqox 



1pR,l,(ri>. 



L.2, 



(41) 



(jg^ ^ e * [#p+ - 0p - +°' ( 6 



while the triplet ones read: 



a 



~j P i( e p+-<t>p-+ e p+- 



-<t><r+)] 



(42) 



(43) 
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Since 



0, we have (e 



i<t>„ 



=/= and by duality triplet correlations. This expression reads: 



/ e ie p „(z) e i6> p _(o)\ 



e~ x '£~ . This property implies that 
all the intrachain superconducting correlations should de- 
cay exponentially along the chains. On the other hand, 
the interchain correlations are reinforced. The physical 
picture is that in this situation fermions of opposite spin 
on each chain are bound together by the exchange inter- 
action. Whether the dominant superconductivity is the 
interchain one or the intrachain one the depends on the 
value of Kff—. Since we have K a - > 1 in our case, the 
dominant superconducting correlations are the interchain 
singlet and the interchain triplet S z — 0. Note that these 
two operators have exactly the same critical exponents. 
In the case of dominant Josephson coupling, the situation 
is reversed. 



e 2ik F , a x e i(9 p+ -(t> p+ ) e icr(9 a+ -<p a +) 



(44) 



This expression again shows power law correlations with 
the critical exponent: 



\{k p+ + k;1 + k 0+ + k-1) 



(45) 



very similar to what has been obtained in Sec. IIVB 21 



B. favoring the Josephson coupling in the coupled 
chain regime 



3. S z = 1 triplet with a 21tF momentum 

It is also straightforward to derive an expression of 
the operator e 2lkF -' jX 4>R.i,<y4'R.2,a associated to the 2fcp 



Let us consider the regime of small interchain hopping 
and assume that now the intra-chain terms are such that 
J|| ~ 4f || . In this case, the dominant term becomes the 
Josephson coupling 5 - and the perturbation term Eq. |36|) 
is replaced by: 



2Ac 



(2ttq) : 



■cosV2(6 p<1 -e p<2 ) 



cos V2(0<t,i ~ 6(7,2) + 77 cos \/2(</><t,i - $0,2) + 7: cos \/2(</><r,i + 0ct,2 + 2wmx) 



(46) 



r 



The treatment parallels the one of the exchange coupling 
in Sec. IV A 11 The Josephson term is rewritten in the 
form: 



2A . 
— — / dxcos26 p . 



cos 26 „ 



- cos 2<p a 



(47) 



The scaling dimensions of the term cos 2# p _ cos 2# cr _ 
is Kp_ + K~_ and scaling dimensions of the term 
cos 2<fi p- cos 20o._ is Kp} + K a - . A similar RG analy- 
sis to the one of Sec. IV All yields for K a - > 1 (0 P _) = 
(mod 7r) and (0 a -) = f (mod 71") and for K a _ > 1 
(0 P -) = (mod 7r) and (</> CT -) = § (mod 71"). This time, 
interchain correlations are decaying exponentially, and 
the dominant superconducting correlations are the intra- 
chain ones. Since K a ^ > 1, (0 a -) = | and the intra- 
chain S z = 1 triplet superconductivity is dominant. For 
Kcr- < 1 the dominant superconducting correlations are 
the intrachain singlet and intrachain triplet S z = 0. 

The results are summarized in the table IIIII where the 
critical exponents r\ a are defined by (O a (x)O a (0)) ~ 
|a;|-'7t«. We note that triplet S z — 1 superconducting 
correlations can never coexist with singlet superconduct- 
ing correlations. 

We have tried to observe these predictions numerically 
but using both the t- J model with small t± , J± and Hub- 
bard model with small t± . Unfortunately, for large J\\/t\\, 
the system with open boundary conditions has edge ef- 
fects due to the proximity of the phase separation which 



occurs generically in the t-J model for large J ft. When 
edge effects are absent and in the Hubbard model, we did 
not find evidences of the proposed predictions, i.e. we 
mostly found cases related to the coupled bands regime 
fixed point. Still, contrary to the isotropic case which has 
been widely studied numerically at zero magnetization, 
a systematic study of the phase diagram at zero magne- 
tization would be necessary before tackling the system 
under magnetic field. This systematic study is out of the 
scope of the present article. 



VI. CONCLUSION 

In conclusion, we studied the case of two fcrmionic 
coupled chains, or ladders, under a magnetic field in- 
ducing a Zeeman effect in the system. The first situa- 
tion we addressed was the free and strong-coupling lim- 
its. We found that large doping-dependent magnetiza- 
tion plateaus occur for "coupled-dimers" systems (with 
large interactions on the rungs) and that pairing is not 
expected in the m = 5 magnetization plateau phase. 
For a system with isotropic couplings, we showed that 
m = 5 plateaus also exist and but pairing survives to 
much higher magnetizations. Furthermore, the compu- 
tation of the one and two-particles gap and bosonization 
interpretation proved that the plateau phase is metallic 
while the system is in a superconducting state below and 
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Dominant interaction 


Dominant superconducting correlations 


Critical exponent 


Exchange coupling (K a ~ > 1) 
Josephson coupling (K a - > 1) 


rung singlet and rung triplet S z = 
leg triplet S z = 1 


i 1 K °+ 

2K p+ 2 

1 1 1 
2K P+ 1 2K a + 



TABLE III: The different dominant superconducting fluctuations with the associated critical exponents. Josephson coupling 
dominates for K p - > 1 and exchange coupling dominates for K p - < 1. 



above this plateau. In addition, we computed the su- 
perconducting upper critical field which is much larger 
than the Pauli limit for a wide range of the parameters 
of the t-J model. Superconducting correlation functions 
precisely agree with the bosonization predictions up to 
the accuracy of our methods. Turning to the coupled- 
chain regime, other interesting unconventional behaviors 
are predicted for the superconducting correlations with, 
for instance, the possibility of having polarized triplet 
correlations under high magnetic fields. However, quick 
numerical investigations study were not able to find good 
parameters providing evidences of such predictions. 

Discussion on experiments - We now briefly discuss 
consequences for experiments. First of all, the possibility 
of measuring irrational magnetization plateau would give 
a direct access to the hole doping S which has recently 
been estimated— to be S ~ 0.1 in the superconducting 
phase of doped ladders. However, the magnetization is 
a global measurement including the contribution of all 
subsystems. For instance, the magnetization curve of 
SCCO has very recently been measured 5 ^ showing that 
the main contribution comes from the chains subsystem. 
On the other hand, measurements such as NMR under 
high magnetic field can provide local information on each 
subsystem. 

The superconducting critical field H C (T) of SCCO has 
been measured under high pressur o 20 ' 21 and displays a 
strong anisotropy and an anomalous H C (T) curvature. 
Furthermore, H C (T) much larger than standard Pauli 
limit are found, suggesting an exceeding of the latest. 
The nature of the superconductivity in SCCO is a long- 
standing debate because of its complex structure and 
the requirement of high pressure experiments. Questions 



such as the nature of the pairing or the dimensionality of 
superconductivity have not reached full agreement yet. 
One possibility is that superconducting fluctuations of 
the RVB type develops in the ladder subsystem leading 
to true superconducting order once these ladders are cou- 
pled through Josephson coupling (T c being controlled by 
this coupling rather than by the magnetic energy scale 
J) . In this case our study suggests that the FFLO mech- 
anism is relevant at the ladder subsystem level and thus 
could be stabilized when coupling the ladders. The other 
possibility would be that the superconducting phase is re- 
ally two-dimensional ) 55 ' 56 allowing other explanations of 
the exceeding of Pauli limit such as triplet pairing— when 
H=0. NMR measurements^ also showed that supercon- 
ductivity survives under rather high field and p— wave 
superconductivity was suggested. Note that our study 
also proposes the possibility for the emergence at high 
fields of pairing channels not present at H = since 
magnetic fluctuations from which pairing originates are 
strongly affected by the magnetic field. Lastly, we men- 
tion that anomalous curvature of H C (T) as found in ex- 
periments were predicted in the mean-field approach of 
FFLO states in quasi-one and two-dimensional d-wave 
superconductors. 4 
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